Dynamical phase quantum thermometer for an ultracold Bose-Einstein Condensate 
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We introduce a primary thermometer which measures the temperature of a Bose-Einstein Con- 
densate ranging from the sub nK regime to the condensate critical temperature. We show, using 
Fisher information, that the precision of our technique improves the state of the art in thermome- 
try. The temperature of the condensate is mapped onto the quantum phase of an atomic dot that 
interacts with the system for short times. We show that the highest precision is achieved when the 
phase is dynamical rather than geometric and when it is detected through Ramsey interferometry. 
Entanglement among several atomic dots can be used to further improve the precision and reach 
the Heisenberg limit. Standard techniques to determine the temperature of a condensate involve an 
indirect estimation through mean particle velocities made after releasing the condensate. In contrast 
to these destructive measurements, our method involves a negligible disturbance of the system. 

Temperature is a crucial concept in quantum physics. Paradigmatic phenomena such as superconductivity and 
Bose-Einstein condensation only occur below a critical temperature. Bose Einstein Condensates (BECs) [1 allow the 
study of quantum effects in systems consisting of up to 10^ atoms by cooling them to regimes in which the individual 
atomic wavefunctions overlap. In this case, the system exhibits quantum behaviour at mesoscopic scales. Very low 
temperatures are also required to observe quantum field theory effects e.g. Unruh-Hawking radiation [2]. 

Quantum field theory predicts that the dynamics of spacetime or the presence of horizons produce quantum particles 
from vacuum fluctuations. For example, in the dynamical Casimir effect, a moving boundary condition gives rise 
to vacuum excitations. This effect was recently demonstrated in superconducting circuits [3 . Currently several 
experimental groups in the field of analogue gravity attempt to demonstrate cosmological particle creation, dynamical 
Casimir effect and Hawking radiation in analogue spacetimes produced in BECs [4 . However, in a BEC a thermal 
background is always present due to unavoidable atomic collisions. In experiments, it is crucial to work at low 
temperatures so that the quantum particles created by effects of emergent spacetimes can be distinguished from this 
thermal noise. Recent experiments to demonstrate the dynamical Casimir effect in a BEC[4 , have not been able 
to ensure that the excitations produced are indeed quantum even at temperatures as low as 200 nK. A method to 
accurately determine temperatures in the nK regime would be of great benefit to these experiments. 

The temperature of a BEC is commonly estimated indirectly by comparing the density profile of the atoms with 
a classical Maxwell-Boltzmann velocity distribution (TJ |5] . This density profile is determined by absorption imaging 
which is a destructive method that involves releasing the condensate. The relative errors in the measurement can be as 
high as ~ 20% [6 . Phase-contrast imaging [7 might be an alternative method to determine the temperature through 
the phase shift of a probe laser beam due to the interaction with the sample. However, at very low temperatures the 
method has low spatial resolution and can only be applied reliably if the condensate is allowed to freely expand for 
a given time [8j. Indeed absorption imaging is the standard technique for measuring temperatures in the nK regime 
[6] and, in general, for probing the condensate [9 . From a more fundamental viewpoint, it is also natural to ask if 
it is appropriate to measure temperatures in a quantum regime using a method which involves a classical velocity 
distribution. 

Recent theoretical and experimental developments in cold atom gases have lead to the experimental demonstration 
of systems consisting of mixtures of two different atomic species or two hyperfine states of the same species p^HT4] . 
The number of particles in each component of the mixture is determined by the interatomic interactions which can 
be modulated experimentally, through a Feshbach resonance [TSHIl]. In particular, atomic quantum dots with only a 
single atom in a tight trap provided by a laser beam can be achieved within a condensate p!6]42Q] . In this paper we 
show how an atomic quantum dot immersed within the condensate can be used to measure the temperature of the 
condensate with high precision and a negligible disturbance of the system. The scheme can be implemented using 
an optical lattice with a single atom per site, where each atom interact at a different time with the condensate [21j. 
This constitutes a primary thermometer which measures the temperatures ranging from the sub nK regime to the 
condensate critical temperature. 
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FIG. 1: a) Sketch of the experimental setup: several atomic quantum dots are embedded within a BEC reservoir. They are 
coupled at different times through a Raman transition to the phononic fluctuations of the BEC. b) Ramsey interferometry 
scheme to measure the dynamical phase and hence measure temperature, c) Ramsey interferometry scheme with entangled 
input states, d) Mach-Zehnder interferometer scheme with entangled input states. A. Input state, B. Interaction switched on, 
C. Final state D. Readout 



Results We start by showing how a two-level system interacting with a quantum field in a thermal state through 
the standard Jaynes-Cummings Hamiltonian can acquire a dynamical phase that depends on the temperature of the 
field. This phase is then read out through Ramsey interferometry. We estimate the precision of this procedure by 
computing the Quantum Fisher Information. Finally, we give details of possible physical implementations of the 
thermometer using atomic quantum dots or an optical lattice coupled to phonons in a BEC (see Fig jl}. In both cases, 
the atoms must interact at different times with the condensate to ensure statistical independence. We show that with 
realistic experimental parameters, our thermometer precision improves on state of the art thermometry after a small 
number of measurements, where in this case, each dot embedded within the condensate or loaded into the optical 
lattice provides one measurement. In the case where the dots are prepared in an entangled state, the precision can 
reach the Heisenberg limit, that is the limit imposed fundamentally by quantum theory. 

We now proceed to introduce our methods and results. We consider the Jaynes-Cummings model that describes 
the interaction of a two-level system of frequency gap Vtd and one mode of a quantum field with frequency Via 

Hjc{t) = frtaO^a + ^^cTz + hg{(j-a^e'^ + a+ae"*^), (1) 

where a, are creation and annihilation field operators, = (1/2) ((Ja^ ^^^y) and (jy are Pauli matrices. The 

coupling strength between the field and atom is given by the detuning is ^ = — and the phase = kx — St 
is a function of time and the position x of the atom. Later on, we will show under which circumstances an atomic 
dot interacting with the phonon field of a BEC can be reduced to this simple interaction. 

The temperature of the field will be determined through the phases acquired by the atom through its evolution 
under the Hamiltonian Q, assuming the adiabatic approximation (see Appendix A). The total phase can be split 
into a geometric and a dynamical part. It has been shown that the geometric phase can be used to gain information 
about the quantum state of a bosonic field. For pure field states, the phase encodes information about the number of 
particles in the field ^22^. In particular, for initial squeezed states, the phase depends also on the squeezing strength 
[23] . More recently, the geometric phase was also employed to estimate the temperature of the field [24]. We will show, 
using Quantum Fisher Information, that the dynamical phase is more convenient to estimate temperatures. We are 
interested in the case where the field is initially in a thermal state pp = "^nPn 1^/) (^/l^ where Pn = {e~-^)^ (1 — e~^). 



3 




FIG. 2: a) Comparison of geometrical (solid, blue) 7^ vs. dynamical 7^ with Qa = 2 tt x 10 Hz, g — 2ti x 0.2 Hz, ^ = 2 tt x 2 Hz 
and c = 5 mm / s. b) Dependence of the dynamical phase Vd with g and ^ at T = 0.5 nK c) Vd (green, solid) and relative 
error 5T /T after 1000 (blue, dashed), 3000 (red, dotted) and 5000 (yellow, dash-dotted) measurements. The parameters are the 
same as in a). 



F = We find, that the geometric phase is 



and the dynamical phase 



- ^^9 e- 

Therefore, the dynamical phase is always larger than its geometric counterpart 

7g<7d- (4) 

In Figj2|i) we plot the dynamical and geometric phases for realistic experimental parameters. One observes that 
for temperatures in the sub nK regime, the geometric phase is negligible in comparison to the dynamical phase. 
Furthermore, we will show that the precision obtained through the dynamical phase in the measurement of the 
temperature is higher than the precision attainable with the geometric one. 

Measurement Since the phases acquired by the system are global and cannot be measured directly, it is necessary 
to setup an interferometric scheme that will allow us to measure relative phases. In particular, the geometric phase, 
used in [24] to estimate temperatures, can only be detected using a Mach-Zehnder interferometer. In this type 
of interferometer the qubit is split into two different trajectories. In this case the qubit can be made to interact 
with a thermal field only in one arm of the interferometer acquiring a relative phase that can be measured at the 
interferometer's output. Now we will use standard techniques of phase estimation in Quantum Metrology to estimate 
the precision of the temperature measurement in this experiment (see Appendix B). 

We consider that in one arm of the Mach-Zehnder interferometer the qubit, initially in state |^), interacts with 
a field in state pr- In the second arm we assume the qubit in state |e) does not interact with the field. Since, up 
to the first order in g/S the state \g) {g\ ® pr — ^Pn ) I? the relative phase at the output will therefore be 
7~ = 7~ + 7~ . Assuming that the dynamical phase can be canceled [25 , we calculate the Fisher information to 
estimate the precision on the temperature measurement through this phase. Since in this case F{pt) = (9t7g)^5 
Cramer- Rao bound in the geometric case employing Mach-Zehnder interferometry, denoted by ST^, yields 

.rs > > (5) 



Here the bound is obtained using a series expansion in g/S. Following the same procedure we find bounds to the 
Fisher information for the case where the temperature is measured using the dynamical phase: 



^/MSirFg'^e^ ^/MAtt F g'^ e^' 
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From Eq. |5]) and (|6| it follows that the error in estimating the temperature using the dynamical phase is smaller 
ST^ < ST^ Furthermore, using the dynamical phase allows us to simplify the setup by using Ramsey rather than 
Mach-Zehnder inter ferometry. Mach-Zehnder atomic interferometry represents a strong experimental challenge while 
Ramsey interferometry does not involve spatially splitting the atomic wave function to undergo different trajectories. 

In Ramsey interferometry the two- level system is prepared in a superposition state, |t) = |^) + |e), and it is allowed 
to interact with the quantum field in a thermal state pr during a time t. After the interaction the probability 
of finding the two- level system in the excited state depends on the dynamical relative phase, F/), picked up by it. 
Under weak adiabatic evolution the state |t) (t| (g) — Zl^n(|^— ) + + (^+1) acquires the relative phase, 

= 1d~^d^ which is very sensitive to the temperature in an ultralow regime such as nK. The Fisher information in 
this case is given by F{pt) = (St7^ —QtId)'^ ~ 4(St7^)^- Therefore, the precision in the temperature measurement 
is increased by a factor of 2 for Ramsey interferometry, i.e. 8T^ = (1/2) JT^. 

In Fig. [2]3 we plot the dependence of F/) on the detuning S and the coupling strength g for the relevant experimental 
regime that we consider below. Notice that the phase is sensitive to g/S and it is essentially determined by this ratio 
and rta- In Fig- ^ we plot the relative error ST^ /T after 1000, 3000 and 5000 measurements. The relative error 
achieved is around 5 %, which improves the 10-20 % of the current state of the art thermometry in the sub nK regime. 

Experimental implementation We now consider a realistic implementation of the above results. We will consider 
a system consisting of a BEC superfluid reservoir in a shallow confining trap interacting with either an array of 
atomic quantum dots or a deep optical lattice loaded with cold atoms I26H28] . Our discussion applies to both, 
dots and atoms in a lattice. However for the sake of clarity in our presentation, we will focus on the quantum dot 
scenarios. In the absence of atomic collisions, the BEC can in principle reach absolute zero temperature. Under 
this assumption, the condensate is commonly described by a classical density function. However, collisions are always 
present and therefore, in the superfluid regime, the condensate is better described by a mean field classical background 
plus quantum fluctuations. The fluctuations, for frequencies lower than the so-called healing length, are given by a 
phononic quantum field operator H(x), which can be expanded in terms of Bogoliubov modes. Now we proceed to 
describe the quantum dots and their interaction with the phononic field. Each dot is created by applying a localised 
steep potential which traps atoms of either a different hyperfine state of the same atomic species as the BEC, or 
alternatively, a different atomic species. We choose a large collisional interaction strength gaa inside each dot well, to 
guarantee that the occupation number in each site is either or 1, giving rise to an array of two level systems. We 
also assume that the wells are deep and separated enough to neglect direct interaction between different sites [26 . 

The phonons and the dots can be coupled through a Raman transition. The effective Hamiltonian obtained under 
the above assumptions can be written as [16l [26]: H = Ha + Hb + Hab- Ha is the free Hamiltonian of an array 
of M dots, which is given by Ha = X^^^i hftd^^l/"^- The frequency ftd is a function of the Raman frequency ft and 
the number of atoms in the condensate Hb is the free Hamiltonian for the phonons Hb = "^Zk ^k^^k^k^ where 
ujk = c |k| and c is the speed of sound. Finally, the interaction Hamiltonian is given by: 



where is the position of the dot in the condensate. The coupling strengths between dot-condensate and condensate- 
condensate atoms are given by gab and ^55, respectively. The frequency is a function of gab^ Qaa^ the Raman frequency 
and the detuning. All the coupling strengths gabidaa and gbb can be tuned with Feshbach resonances. In particular, 
we will consider that the strengths are tuned such that 5' = {). In this case the system is described by the spin- 
boson Hamiltonian [16] , in which the structure of the reservoir is characterised by the spectral density. The ubiquity 
of physical systems that can be described by the spin-boson Hamiltonian directly illustrates that the thermometry 
technique we describe can be applied to a range of different physical systems. 

Now, with suitable boundary conditions for the condensate trap -e. g. hard-wall or box potentials [29l [30]-, the 
energy separation of the phononic modes can be large enough to ensure that each two-level system is only effectively 
coupled to the mode with closest frequency to Vtd- We denote this frequency by and apply the Rotating wave 
Approximation (RWA). Assuming that each dot interacts with the condensate at different times, the Hamiltonian for 
each dot is given by Eq. Q, where now the dot-field coupling strength g is given by: 



M 



M 




(7) 




(8) 



where V is the condensate volume. The RWA holds for weak couplings and small detunings {g^d} « e where 
e = Qa ^d- Taking typical values L = 500 jam and c = 5 mm /s and assuming \f the frequency of the phonons 
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is given by = 27r x 10 Hz. With experimental data for the Feshbach resonances p!3HT5] the values of g can be 
tuned within a broad range from 27r x 0.1 Hz to 27r x 10 Hz. Taking g around 0.1 Hz and a detuning of 27r x 1 — 2 Hz 
we are consistent with both RWA and the adiabatic approximation. This is the range of parameters explored in the 
plots. 

Discussion The description of the experiment is the following. The atomic dots are prepared in a separable 
superposition of their internal levels, 0^=1... m |ti)- They are then coupled to the phonons through Raman transitions 
and allowed to interact with the BEC for some time, which in our plots is of the order of O.ls. The times can be made 
significantly shorter by relaxing the cyclicity condition, which we imposed in order to compare the dynamical phase 
analysis with the geometrical phase analysis. However, it is no longer necessary in this setting. The probabilities of 
excitation are then measured using standard imaging techniques j^. In order to make sure that the measurements 
are completely independent, each dot must interact at different times with the condensate. Moreover, the separation 
in time between interactions should be large enough to rule out temporal correlations between the occupation numbers 
of the considered mode. These temporal correlations drop to after a few ms for typical parameters [31 .Another 
way to ensure that different measurements are not correlated is to couple each dot to a different mode of the field. In 
this case, the interactions can be simultaneous, since each dot would probe a different thermal distribution with the 
same temperature. The only limitation to the number of probes is that the spatial separation must be larger than 
the healing length of the condensate in order to ensure no direct interaction between them [26 . Moreover, indirect 
interaction due to phonon exchange is proportional to g'^/S^ and therefore is negligible for times t = 1/S and g/S « 1. 
Taking for instance a realistic separation between the dots of 300 nm and a condensate length of 500 /im, we can load 
1500 impurities in an optical lattice. Since the lifetime of the condensate can be larger than 100 s it seems feasible to 
perform at least a few measurements on each dot. Therefore, as can be seen in Fig. ^) the relative error is well below 
the 20 % state of the art ultracold thermometry [6 , and ultrahigh precisions below 1 % can be achieved. Moreover, 
as explained above, by preparing several dots in an entangled input state [18» i32j the precision would be significantly 
improved reaching the Heisenberg limit. 

The Unruh temperature can in principle be measured in our setting by accelerating the dots. The idea of demon- 
strating the Unruh effect in BECs was suggested with a different setup in [33 . Moreover, the use of a geometric 
quantum phase as an estimator of the Unruh temperature was proposed using a challenging experimental setup 
involving atom Mach-Zehnder interferometry [34]. We are currently working on implementing our technique to mea- 
sure the Unruh effect in a BEC with higher precision. If a is the acceleration, the Unruh temperature is given by 
Tu = h a/ {27r cks) which in our setup give rise to Tu = 2.4 nK for an acceleration of g -where g is the acceleration of 
gravity on the surface of the Earth. Accelerations of several g were achieved in recent experiments [21 . Therefore, 
with our method the detection of the Unruh temperature of the phononic BEC bath is clearly within reach of current 
technology. 

We have shown that the quantum dynamical phase acquired by a two-level system interacting with a quantum 
field in a thermal state under the Jaynes-Cummnigs Hamiltonian can be used to measure the temperature of the 
field. Since the Jaynes-Cummings model successfully describes a number of different quantum systems, our ideas 
can be applied to measure the temperature in a great variety of physical setups, such as cavity QED, ion traps or 
circuit QED in a weak coupling regime. As a particularly interesting implementation we choose to present in this 
paper a cold atom setup. In particular, an atomic quantum dot is coupled to the phononic quantum fiuctuations of 
a BEC. We compute the dynamical phase acquired by several independent dots in a Ramsey interferometry scheme 
and estimate the precision by means of the Quantum Fisher Information. We show that the phase is sensitive to 
temperatures in the sub nK regime and that the precision improves on state of the art thermometry after a number of 
measurements that can be achieved with current technology. Moreover, the precision can be improved by preparing 
entangled input states. In summary, we have introduced a non-invasive, accurate and quantum way of measuring 
ultralow temperatures in a Bose-Einstein Condensate. In addition to improving on state-of the art thermometry, 
these techniques provide an experimentally feasible method of measuring very-low temperatures arising in analogue 
gravity scenarios, paving the way for the detection of analogue Unruh and Hawking temperatures in BECs. 

Appendices 

Appendix A: Dynamical and geometrical phases The time evolution of the system can be found by writing 
Hjc{t) = U{e)H%U{0)^ where U{0) = e^^^^'^ and H% = Hjc{x = 0,t = 0). The eigenstates of are given 

by 1^—) = cos(^) |(n + 1)^) — sin(^) |ne) and |n+) = sin(^) |(n + 1)^) + cos(^) |ne) where \ng) and \ne) are 

the eigenstates of the free Hamiltonian (corresponding to ^ = 0) and an = arctan( ^^ ^^^"'^ ). In this notation |n) are 
field number states and \g) and |e) are the ground and excited states of the qubit. We require that the qubit only 
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acquires phases during its evolution. Therefore, to ensure that there are no transitions between the U{0) \n±) states, 
we make use of the adiabatic approximation [35 ^^^5^^^ t << 1. In the Jaynes-Cummings model transitions 

are only possible between the eigenstates with same n, therefore the adiabatic condition yields ^-^^^^j^^^—^t « 
For simplicity, we assume that the parameter undergoes a cycle, therefore the interaction time is t = 2t:/5. The 
dynamical phase acquired by the state |n±) under this condition is given by 

= -y(fi.(..-i)±\/i2 + 4g2ll). (9) 

The geometric phase defined by [22j: i7g„ = < U^dgU >„± d9, yields 'Jq^ = 27r(n — cos^(^)) and = 
27r(n-sin2(^|i)). We are interested in the case where the field is initially in a thermal state pF = Xln^^ 1^/) (^/l' 
where Pn — (e~^)^ (1 — e~^), F = In the mixed case, the geometric phase is given by 7^ = ^^s(Xln^^ e^^^c^)) 

[36l [37] , while in the dynamical case, the phase is obtained through a simpler expression 7^ = PnlDw 
find, using a series expansion on g/S^ that the geometric phase is given by Equation ([2| and the dynamical phase by 
Equation ([3| 

Appendix B: Quantum Fisher Information The quantum Cramer- Rao bound [38^ states that the error Sx 
in estimating a parameter x with M measurements on a state px that depends on x, is bounded by: 

5x > , ^ (10) 

- ^yuFipxY 

F{px) being the Quantum Fisher Information of the state. This inequality assumes that input states are separable, 
however, it has been shown that entangled states improve the shot-noise scaling reaching the Heisenberg limit 

given by a 1/M scaling where M is the number of entangled particles. F{px) can be written as [39] : 



(m| dxpx \n) 

Pm + Pn 

where pm,n are the matrix elements of px in an orthogonal basis in which px is diagonal: px = Pn 1^) (^1- The 



F(p.)=2^ l^"^l"-f (11) 

Pm ~r Pn 



terms in which both and pn are zero in the sum in Eq. (11) must be excluded. 
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